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Abstract 

We introduce a notion of universality classes for the Gregory-Laflamme instability and deter- 
mine, in the supergravity approximation, the stability of a variety of solutions, including the 
non-extremal D3-brane, M2-brane, and M5-brane. These three non-dilatonic branes cross 
over from instability to stability at a certain non-extremal mass. Numerical analysis suggests 
that the wavelength of the shortest unstable mode diverges as one approaches the cross-over 
point from above, with a simple critical exponent which is the same in all three cases. 
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1 Introduction 

The Gregory-Laflamme instability |2] has been conjectured to arise precisely when local 
thermodynamic instabilities exist, for horizons with infinite extent and a translational sym- 
metry Pill]- The idea is that, barring finite volume effects or some unexpected consequences 
of a curved horizon, the horizon will tend to become lumpy through real-time dynamics pre- 
cisely if it can gain entropy by doing so. Based on this line of thinking, one might poetically 
ascribe the dynamical stability of the Schwarzschild black hole in four dimensions merely to 
the fact that the horizon is smaller in total size than the shortest wavelength instability. 

A fairly thorough analysis in exploiting a connection between time-independent, 
finite-wavelength perturbations in Lorentzian signature and negative modes in Euclidean 
signature, gives considerable confidence that the conjecture of O EJ is correct, as well as 
making clearer the reasons why thermodynamic and dynamical stability should be connected. 
The analysis of j3] does leave several points open, including the following: 

1. For the black brane solutions familiar from string theory, just when is stability lost? 

2. What happens near the threshold of stability? 

3. Once an instability develops, what does the system evolve into? 

Question 3) has been the subject of recent scrutiny for the case of uncharged black branes 
(see for instance [HJ |Hl El EI])- Question 1) is straightforward if one grants the analysis of 
jS], since one needs only compute the specific heat as a function of non-extremality. We will 
address the question in section |21 from this thermodynamic point of view, and in section 0] by 
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directly searching for time-independent, finite- wavelength perturbations, thus providing an 
independent check on the validity of the analysis of We will also make a start on question 
2), based on our numerics. In section |3] we will introduce a notion of universality classes 
for black brane horizons, based on solutions described in [TT]. For the questions we will 
address, a single characteristic exponent defines a universality class: this exponent describes 
how horizon entropy scales with temperature. 

Our notion of universality classes is of some intuitive use in understanding the stability 
properties of various charged brane solutions. In particular, the well-known non-extremal Dp- 
branes of type II string theory interpolate between one universality class far from extremality 
(where the entropy scales as a negative power of temperature, and Gregory-Laflamme in- 
stabilities exist) to a different universality class near extremality (where entropy scales as a 
positive power of temperature, and Gregory-Laflamme instabilities probably do not exist). 

Other recent work on the Gregory-Laflamme instability includes [T2*l FEU IT^j . In particu- 
lar, ^2] explores the stability of charged brane solutions, using methods somewhat different 
from ours, but with conclusions that overlap and agree with those in section 0] 



2 Thermodynamic considerations 

The Dp-brane solutions to type II string theory have the following string frame metric and 
dilaton: 

ds 2 = G tt dt 2 + G rr dr 2 + G ee dtt 2 n+1 + Gyydtfdtf 

h d t> + ^Z rfr 2 + Jfr 2 dn 2 n+1 + -^=dy l dy l (1) 



e 24> _ y(n-4)/2 



where 



rn sinh 2 a 



h=l-± / = l+ ' 0&lim (2) 



The string metric G ^ is related to the Einstein metric g^ v by = e~^/ 2 G^ v . In equa- 
tions and (J2J) , n — 7 — p. We will assume throughout that the metric ((TJ is being used to 
describe a large number N p of coincident Dp-branes, so that the supergravity approximation 
is reliable (except perhaps near the horizon at extremality, where some of the solutions have 
a null singularity). 

The thermodynamic properties of Dp-branes are well-known, so we will present here only 
a brief summary. The original literature, including ^1^21 El, can be consulted for further 
details. The entropy of the Dp-brane solution (JTJ is 

„ _ ^Einstein _ ^string _ 4irQ n+1 n+1 , , . 

b ~ AG N ~AG N e^~ 2k 2 Vr ° C ° Sha ' [6) 

where ^Einstein is the horizon area measured using the Einstein frame metric, A str - mg is the 
horizon area in the string frame metric, V is the coordinate volume in the y l directions, and 
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f2 n+ i is the volume of the sphere S n+1 : 

m + l 

n m = Vo\S m = 2 - 2 . (4) 

r 

If one compactified the on a torus whose volume at infinity was V, then (jSJ) would be 
the total entropy of the resulting horizon. The ADM mass (most easily calculated using the 
Einstein frame metric) is 

M = ^ r «(— + 2 C08h2a )' (5) 
The temperature can be extracted from the surface gravity: 



T _ ^_ _ J_ dg t t/dr 



2tt 2tt 2y/g tt g rr 



(6) 



t=tq 



47rro cosh a 



One can substitute G i4 and G rr for and g rr in (0), provided the dilaton is non-singular 
at the horizon (which would just mean that the horizon is a null singularity, and that only 
happens for extremal Dp-branes). 

It is straightforward to check that T = DM/ OS. 

The Dp-brane charge is quantized, and can be expressed in terms of the following con- 
straint: 



2k 2 N p J2K 2 (2iry- 2 P(a>) 3 -P (2tt)" 

R n = |rg sinh 2a = —— N p r p = ^ K —— = v^ 7 ^iV P , ( 7 ) 

ni2„ +1 ni2 n+ i nll n+1 



where N p is an integer and t p is the tension of a single extremal Dp-brane, and we have used 
the relation 

16ttG n = 2k 2 = (2tt) W) 4 , (8) 

valid in uncompactified type II string theory. The formula (J7J), taken from |16j . can be 
checked against the standard formula for D-brane tension: 

2K 2 T 2 p = (2n) 7 ~ 2p {a'f~ p . (9) 

Holding N p fixed, the behavior of the temperature in different limits is well known. In 
all cases, T — > in the small a limit, since this is the large mass limit where the charge is 
negligible. For p = 1,2,3, and 4, T — ► also in the extremal limit, whereas for p = 5, T 
approaches a constant, and for p = 6, T diverges. There are no other remarkable qualitative 
features of the dependence of T on mass other than a consequence of what we have already 
said: for p < 5, there is a special mass M* where T reaches a global maximum. This is the 
border of thermodynamic stability: for masses less than M*, the specific heat is positive, and 
for masses greater than M* it is negative. Holding N p fixed, one finds that (jUJ) is extremized 
when 

sinh 2 a = sinh 2 a* = , (10) 

n — 2 
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or alternatively when 

M M* n 2 -2 

N p TpV N p t p V n^^l ^ ' 

(note that N p r p V is the extremal mass). Another characterization of the same condition is 

^H^r*- <i2) 

where R is the length scale introduced in (JZJ). 
The M2 and M5-brane metrics are given by 
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ds li2 = /m2 /3 (-hhndt 2 + dx 2 2 ) + + r z dVL z 7 



ds M5 = /m5 /3 (~ h M<odt 2 + dx 5 2 ) + f ^— + r 2 dfi 



(13) 



where / and h are given by (J2j) with n = 6 for the M2-brane and n = 3 for the M5-brane. 
Indeed, the thermodynamics for the M2-brane and the Dl-brane are identical. This is not 
accidental: wrapping an M2-brane around a circle gives a fundamental string of type IIA, 
which is equivalent through some further dualities to the Dl-brane (and, more to the point, 
has the same Einstein metric up to a slightly different identification of parameters). For 
similar reasons, the M5-brane and D4-brane have the same thermodynamics. In particular, 
the condition for the maximum Hawking temperatures for the M2-brane and the M5-brane 
are given precisely by (JTHJ), (|TT|) . and fT2"J) . with n = 6 and 3, respectively. 



3 Universality classes and the general problem 

A fairly broad statement of the question of p-brane stability in the supergravity approxima- 
tion can be summed up as follows: starting with some p-brane solution in supergravity, one 
wishes to consider perturbations that break some part of the translational invariance and 
determine whether any of them grow with time. Such considerations can often be reduced to 
a problem in 2 + 1 dimensions, and in the following paragraphs we shall give more particulars 
about how this comes about. 

One starts with a two-derivative action in D > 2 spacetime dimensions, containing Ein- 
stein gravity and perhaps also form fields and scalars. The brane in question is some solution 
with a horizon and a weakly curved region far from it (typically asymptotically flat space), 
where the whole solution admits as symmetries the Euclidean group of translations and ro- 
tations in p dimensions. 1 There should also be some Killing vector which is timelike in the 

Smaller symmetry groups can arise. In "smeared" brane solutions, for instance, composed of a density 
of p-branes distributed in several orthogonal directions, the symmetry group would be translations in all 
dimensions plus a product of two rotation groups for the directions parallel to and perpendicular to the 
individual p-branes. 
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weakly curved region and null at the horizon, which can be used to construct a notion of 
energy; and there may be assorted spacelike Killing vectors as well, indicating rotational 
symmetries. The brane is of spatial co-dimension D — 1 — p, and far from the brane one 
generally expects an approximate SO(D — 1—p) rotational symmetry group. The brane may 
have various conserved charges and angular momenta (the latter corresponding to generators 
of SO(D — 1—p)). In this rather general context, we would like to ascertain the stability 
properties of the solution against small perturbations (determined by thermodynamic sta- 
bility, according to jHlSE]), perhaps also the dispersion relation for such perturbations, and 
finally the end state of evolution along an unstable direction — though this last is surely a 
much harder problem in general than the other two. 

Although we believe that the notion of universality classes we introduce below will have 
some applicability to the general case (including for example smeared brane solutions), let 
us focus on the case where the SO(D — 1—p) symmetry is exact. Then perturbations can 
be labeled uniquely by their wave-number k along the brane and their angular momentum 
quantum numbers under SO(D — 1—p). It seems reasonable to assume that the s-wave 
perturbations become unstable first, provided they exist as dynamical perturbations (as 
opposed to being perturbations which are pure gauge). The perturbations in question thus 
preserve a large subgroup of the symmetry of the original solution: SO(D — 1—p) times 
that part of the Euclidean group that preserves k. This makes it particularly natural to 
consider a Kaluza-Klein "reduction" where one integrates the action over an S D ~ 2 ~ P orbit 
of SO(D — 1— p), 2 and drops the integration over an R p_1 perpendicular to k. The result is 
a 2+1-dimensional action, which in a fairly broad set of circumstances can be brought into 
the form 

S = jd z x^(R- -G^d^WcjP - V($f) . (14) 

Einstein gravity, form fields, and scalars in D dimensions reduce down to the type of action 
written in (JHJ), plus perhaps abelian and non-abelian gauge fields. The abelian gauge fields 
can be dualized to scalars, unless there are Chern-Simons terms in the action. Non-abelian 
gauge fields in general cannot be dualized, but examples in the literature often have fields 
excited corresponding to an abelian subgroup. Thus (|T1J) . while obviously not completely 
general, does cover a broad range of cases. Furthermore, it has been observed [IS] that 
V(4>) is commonly a sum of exponentials of canonically normalized scalars: for example, the 
integrated curvature of S D ~ 2 ~ P scales as a power of its radius, but the canonically normalized 
scalar measuring the size of S D ~ 2 ~ P is some multiple of the logarithm of the radius. 

In section ETT1 we will consider properties of solutions to (|14j). in part recapitulating argu- 

2 An interesting aside is that in circumstances involving angular momentum, where obviously SO(D~l—p) 
is not preserved even by the original solution, Kaluza-Klein reduction is still often useful, with the somewhat 
mysterious "consistent truncation" ansatze providing an exact lower-dimensional account of both the solution 
and some low partial wave perturbations J3J 0]. We suspect that the coincidences arising in consistent 
truncation that matter for analyses of horizon perturbations arise from bosonic symmetries rather than 
supersymmetry. 
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merits of jTH]. In section l3~21 we briefly remark on the reason to search for static perturbations. 
Then in section 13.31 we will make an explicit numerical study of stability of a special case of 
the solutions considered in section 13.11 



3.1 Solutions to gravity plus scalars 

Let us now study the slightly more general situation where 

S = j d D x^g(R- l -G l3 {$)d^d^ - , (15) 

where D = d + 1 is unrelated to D. Motivated by the previous discussion, we will be most 
interested in the case D = 3. Solutions to the equations of motion of this action were studied 
in ^S], of the form 

ds 2 = a{r) 2 {-dt 2 + dx 2 ) +dr 2 , $ = ${r) , (16) 

with a(r) a monotonically increasing function of r, running from to oo. The large a region 
in these solutions is weakly curved (in fact, the main interest in [T^j was asymptotically 
anti-de Sitter space), and the small a region is in most cases singular. These solutions have 
no horizons (or at best degenerate or singular horizons), and their symmetry under boosts 
amounts to the statement that they are extremal. All such solutions can be derived using a 
first order formalism (see for example [T9~]I2UI 12T] ). which is a direct analog of the Hamilton- 
Jacobi method of generating FRW cosmologies [22] : if one starts with an action of the same 
form as (|14|) . but in D dimensions, and finds W((j)) such that 

V{<P) = 2° WW ~ 4D~2 W ' (17) 



then a solution to 



*Z-W=fF® (18) 

will also solve the equations of motion. 

Non-extremal generalizations of (jlfij) have the form 

dv^ -» -* 

ds 2 = a{r) 2 {-h{r)dt 2 + dx 2 ) + j-- , = 0(r) . (19) 

ri[r) 

We assume in general that h(r) — > 1 in the region where a(r) — > oo, and h(r) < 1 everywhere. 
If h(rn) = and r# is the largest value of r for which this is true, then r = rn is an event 
horizon with respect to the region where a(r) — > oo (which is also at large r). Provided 
h'(rfi) 7^ 0, this is a non-degenerate horizon with finite Hawking temperature. 

Solutions of the form (|19|) involve a number of constants of integration. The counting of 
them has been explained in [IB], and we will recapitulate briefly here. The relevant equations 
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of motion are all ordinary differential equations in r once we assume the form ()19|). They 
comprise the second order equations for the scalars, of which we suppose there are n; the G rr 
component of Einstein's equation, which is a first order equation; and one further Einstein 
equation (for instance, the combination G t t — G x x of Einstein's equations), which is second 
order. There are 2n + 3 integration constants altogether. One amounts to an additive shift 
in the radial variable r, which is just a coordinate transformation. Another is fixed by 
the requirement that h(r) — > 1 in the region where a(r) is large. The presence of a non- 
degenerate horizon fixes n more, as explained in [T^] (one way to understand the presence of 
these n horizon constraints is that in the Wick-rotated Euclidean solution, where the horizon 
becomes a point, the scalars must be smooth everywhere, so their radial derivatives must 
vanish at this point). This leaves n + 1 parameters that specify the solution. One of these 
is the temperature of the horizon, and the other n pertain to the asymptotics of the scalars 
in the large a region. Of these n parameters, some but not all may be fixed by demanding 
regularity in the large a region. In an AdSd+i/CFTd context, where the large a region is 
asymptotically anti de-Sitter, the n parameters correspond to the coefficients of e~^ d ~^ r l L 
in the expansion of the scalars around their constant limiting values as r — > oo. These 
coefficients amount on the CFT side to the mass parameters of gauge singlet operators added 
to the lagrangian. If some of the operators in question are irrelevant, then the corresponding 
scalars have positive m 2 , and the larger of the two solutions to the linearized scalar equations 
of motion blows up at large r: this is a case where parameters are fixed to zero by demanding 
regularity in the large a region. In the case where the dual operators are relevant, and the 
corresponding scalars have negative m 2 , then the coefficient of e -( d ~ A i) r / L [ s truly a variable 
parameter. In the case of asymptotically anti-de Sitter solutions, an additional scaling 
symmetry can be used to eliminate one additional parameter: this amounts to applying a 
scale transformation on the CFT side to set, say, the temperature equal to 1. 

Let us now consider the situation where the asymptotics of the scalars is entirely fixed in 
the large a region, leaving us with one free parameter. One can show (see [T^j for the case 
where d = 4) that 



and then B is the free parameter. (The identity (|20|) holds even though a(r) may in general 
change when B changes. It is a convenient parametrization because the influence of B on 
a{r) is small in the large a region). In an AdSd+i/CFTj context, this would correspond 
to specifying a deformation of the gauge theory lagrangian by relevant operators and then 
varying either the temperature or the energy density. Evidently, in the limit B — > 0, one 
recovers a solution of the form ()16|) . Suppose that in this limiting solution, V(<f)(r)) decreases 
monotonically as a(r) decreases, going to oo as a — > 0. 3 It can then be argued (though not 
wholly rigorously) [TB] that 

3 It is possible that this condition could be relaxed somewhat: we must certainly have V{<j>{r)) less than 
its limiting value at large a, and probably the monotonicity property is also needed for a less than a certain 
upper bound. 




(20) 
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• Solutions with regular horizons exist for all B > 0. 

• If V(4>) is the sum of exponentials of canonically normalized scalars, then in the B = 
solution, the scalars have an asymptotic direction as a(r) — > 0: that is, <f)(r) ~ 
0o + vip(r) for some constants <fro and v, which we may normalize so that tp is itself 
canonically normalized. Then V(4>(r)) ~ e 2llp in the region of small a. 

• In near-extremal solutions, with sufficiently small B, the solution at small a is well- 
approximated by a Chamblin-Reall solution for the exponent 7, which we explain below 
in section 13.31 The whole solution can be well-approximated by the original extremal 
solution patched onto a Chamblin-Reall solution. These approximations become pro- 
gressively better as B — > 0. 

• The entropy and temperature of near-extremal solutions are related by a power law: 
S oc VT a , where 

a= l-2 7 *(D-2) - (21) 

These features suggest the notion of universality class that we are going to use. Suppose 
that an extremal solution to an action with canonically normalized scalars has a "scaling 
region" where 0(r) pa </> + vip(r) and V ~ e 2w for many decades of variation of a(r). Then 
non-extremal solutions with the horizon well within the scaling region should again be well- 
approximated by the original extremal solution patched onto a Chamblin-Reall solution near 
the horizon. The approximation should become good in the limit where the scaling region 
in the extremal solution extends far above and below the value of a where the horizon is 
located in the non-extremal solution. The same power law, S oc VT a , should pertain. 

With the idea of universality classes in hand, we can address stability of the horizon in 
a simple way. In section 13.31 we will study explicitly perturbations around Chamblin-Reall 
solutions, involving the metric and the "active" scalar tp. We will find perturbations which 
are normalizable both at the horizon and far from it. As the full solution is well approximated 
by matching the extremal solution onto a Chamblin-Reall solution, perturbations should be 
well approximated by similarly matched solutions — only, because the perturbations we are 
considering are normalizable within the Chamblin-Reall region, they approach zero near the 
matching region, and the perturbations of fields in the extremal region are very small. It 
does not seem plausible that there are normalizable perturbations which are large in the 
extremal region and small as one enters the Chamblin-Reall region near the horizon. To sum 
up, we believe the stability properties of the whole solution are determined by the stability 
properties of the near-horizon Chamblin-Reall solution, which also determines the scaling of 
entropy with temperature. 

It may seem that we are focusing rather narrowly on a rather special class of pertur- 
bations: not only invariant under the SO(D — 1 — p) that we started with in carrying out 
the Kaluza-Klein reduction, but also involving only the metric and the active scalar near 
the horizon. This amounts to focusing on a type of fluctuations which we might describe 
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as adiabatic, since the fields are varying locally in the same proportions that they would do 
globally if we simply changed the temperature. But this is precisely the mode of instability 
that the thermodynamic arguments of [21 IH El suggest. Instabilities in other modes are 
conceivable, and one could even study them by considering scalar fluctuations "transverse" 
to the solution. But we would find it very surprising if such fluctuations gave rise to normal- 
izable instabilities when the fluctuations we consider do not. To put it another way, for that 
to happen, the logic of jH] must fail. Our rather detailed predications about the nature of the 
instability amount in a sense to an elaboration of the connection between thermodynamic 
and dynamical instabilities. 

3.2 Static perturbations 

In the rest of this paper, we will investigate Gregory-Laflamme type instabilities of various 
black branes. For all the unstable black branes Gregory and Laflamme analyzed |2] , they 
found at linear order in perturbation theory that there is a static perturbation, and that 
instabilities occur at longer wavelengths than this static perturbation. 

As remarked in the introduction, Reall [S] in his approach to the proof of the conjecture 
°f 13 El investigated the role of such time-independent perturbations, and argued that they 
are of central importance to the relation between thermodynamic and dynamical instability. 
In particular he related such perturbations of a black-brane solution to the negative modes of 
the Lichnerowicz operator on the Euclidean black hole background found by compactifying 
the black p-brane on a p-torus and doing a Wick rotation. The existence of such negative 
modes is then related to thermodynamic stability. (Negative modes of Euclidean black holes 
have also been analyzed in the context of finite temperature stability of flat spacetime by 
Gross, Perry and Yaffe |23j.) 

Thus when looking for instabilities of black branes, instead of considering instabilities 
with arbitrary growth rates, we will restrict our attention to static threshold perturbations. 
With the relation between thermodynamic and dynamical stabilities in mind, we expect the 
dispersion curve of a Gregory-Laflamme type instability to end at such a threshold point at 
the small wavelength limit. 

3.3 The Chamblin-Reall solutions 

As argued in section 13.11 close to the horizon, black brane solutions fT] for gravity coupled 
to a scalar with a potential of the form V = V^e 2 " 1 ^ provide good approximations to a 
broad range of supergravity solutions, including ones that are non-extremal versions of those 
that represent RG-flows in AdS/CFT. 

The action for the gravity plus scalar system is 



,5' 




(22) 



V(<j>) = V e*>+. 
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We will restrict our attention to negative Vq, since that is the case that admits black hole 
solutions. The black brane solutions are given by 



dsl = a>l(r)(—h (r)dt 2 + dx D _ 2 ) + 



dr 2 



(23) 




where 



"»w = i - (?) 



i 



(24) 



a (r) = r 2 ^ 2 ( D - 2 ) 



0o(O = log(7 2 cr) , 
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and c is given by the positive root of the equation 



V4 = C(^ 



ID 



1 



) 



(25) 



4D-2 



In the notation of [TT] (differing slightly from ours), these solutions correspond to the "type 
II" case, with b 2 < D — 1, M > 0, and Vq < 0. Note that in the gauge presented above, 
the metric is independent of Vq: it only depends on 7, whereas the scalar profile depends on 
both 7 and Vq. 

With respect to the general discussion of non-extremal horizons in section l3~Tl the solution 
(124(1 is rather special, in that the functions a(r) and 0(r) are entirely independent of the 
parameter tq that determines the location of the horizon. (The equivalent parameter in 
section I3~T1 is B). In other words, one can get at the solution ()24j) by first obtaining the 
extremal solution and then using (|20|) without changing a(r). In general, ^ 1 "back- 
reacts" on the solution, changing a(r) and <p(r). This property of (|2"4*|) also hints at the 
origin of the relation (|23|): One plugs W = ce 1 ^ in (fT7|) to obtain V, just as one would do 
when looking for extremal solutions. 

As noted in ^T], for a discrete set of 7's it is possible to arrive at the solutions (|24[). 
(1231) by a dimensional reduction of spacetimes of the form Schwarzchild q+ 2 x i? D ~ 2 on the 
g-sphere of Schwar zchild q+ 2- For such cases, V and 7 are given in terms of q as 



At least for this discrete set of solutions, we certainly expect to find dynamical instabilities: 
they correspond to the instabilities of an uncharged black string in q+3 spacetime dimensions. 
The 7's that can be obtained by a dimensional reduction on a g-sphere, where q runs from 
1 to 00, fall in the range 



q(q- 1) 




(26) 




(27) 
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(The case q = 1, which gives the upper limit on 7, is degenerate: V = for that case.) 
It turns out that (J27|) is precisely the range for which the .D-dimensional black brane has a 
negative specific heat. The conjecture of [Sill] then implies that one should have dynamical 
instabilities for the 7's precisely in this range, including the ones that can't be obtained by 
dimensional reduction. In the next section, we will give numerical evidence in favor of this 
for the case D = 3. (The same conclusion for other D then follows from straightforward 
Kaluza-Klein reduction) . 



3.4 Numerical study of the stability of Chamblin-Reall solutions 

Choosing D = 3, we searched numerically for threshold dynamical instabilities of the 
Chamblin-Reall solutions by restricting attention to perturbations of the form 



ho(r) J (28) 



(r) = 4>o(r) + A cos kx<j>\(r) 



where 



(29) 



A(r,x) = XcoskxAi(r) 
B(r, x) = A cos kxBi{r) , 

and working to linear order in A. It is possible to show that this ansatz for the perturbations 
is consistent with the equations of motion. 
The Einstein equations in D = 3 read 

K = TS - gZTZ . (30) 

Four of these equations are nontrivial, namely the ones involving i?*, R x x and R T X . Using 
the R r x equation, one can solve algebraically for B\ in terms of A%, A[, and <pi. The R\ 
equation and the scalar equation of motion involve B\, but not its derivatives. Plugging in 
Bi as found from the R r x equation into the R\ equation and the scalar equation of motion, 
we end up with the following two equations involving A\ and 0i: 

0= (-^k 2 r 2 ) +r^ 2 - 2 7 2 r^ 2 + 7 4 r^ 2 (l + k 2 r)) A x (r) + 7 (-1 + 7 2 )' r r " ^(r) 

+ I 7 2 (_3 7 4 r 2 + r 7" 2 (_ 4r + r 7- 2 ) _ 2 7 2 r 7 " 2 (-5r + 2r^ 2 )) A[{r) 

— 7 4 r (r — r 7 (y 2 r — r 7 A"(r) 
= 2 (-1 + 7 2 ) 2 r 7 " 2 Ai(r) + 7 (2 r 7 " 2 - 7 2 r 7 " 2 (4 + k 2 r) + 7 4 (k 2 r 2 + 2 r 7 " 2 )) <f> x {r) 

- 7 2 (-r + r 7 ~ 2 ) (-3r 7 " 2 + 7 2 (-r + 4r 7 ~ 2 )) A[{r) + 7 (- (V r) + r 7 " 2 ) 2 0' x (r) 

+ 7 3 r (— r + r 7 ) ( — (l 2 rj + r" 1 J ^'{ix) ■ 

(31) 
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Figure 1: Stability of Chamblin-Reall solutions. The dots represent threshold instabilities of 
the Chamblin-Reall solutions. The vertical line at 7* = 1/ v2 indicates the thermodynamic 
stability limit. 

(It is possible to do a consistency check by arriving at the same equations using the R* and 
R r r equations). 

We used the linearity of the equations to fix the value of <pi at the horizon. Then, 
using a a shooting algorithm, we looked for values of k and Ai(r ) that would give solutions 
that are well-behaved at the horizon and infinity, for a given 7. The results for k as a 
function of 7 are shown in figure H As can be seen from the figure, the wavelengths of the 
threshold instabilities diverge as one approaches the thermodynamic stability limit given by 
7 = 7* = (the lower limit in ()27j) ). 7 = 1, which corresponds to the upper limit in (J27j) . 
is where the causal structure of the Chamblin-Reall solutions changes: in the notation of 
Chamblin and Reall , going from 7 < 1 to 7 > 1 corresponds to going from the b 2 < D — 1 
case to the b 2 > D — 1 case of the type II solutions. 

For 7 < 1/ v^2, we found no static perturbation. This plus the thermodynamic stability 
of the solutions strongly suggests that they are stable. 

In order to address the question of what happens near the threshold of instability, we 
attempted to fit the behavior of A; as a function of 7 to a power law near 7 = 7*. The fit 
was considerably better if the threshold of stability was shifted to 7^ um ~ 1.0037*: then we 
found k ~ — Inum- The fit to a power law was still less than spectacular. We suspect 
that numerical error contributed both to this and to the discrepancy between 7* and ^f num - 

4 Crossing the threshold of stability: non-dilatonic 
branes 

Next, we would like to investigate the dynamical stability of D3, M2 and M5-branes. Fol- 
lowing the philosophy introduced in the beginning of section El we will do a Kaluza-Klein 
reduction of the metrics (P) and (|13j) on all but one of the spatial worldvolume directions 
and S n+1 , to end up with a 3 dimensional background. We will be working with s-wave 
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perturbations of this background, keeping only the two scalars that describe the sizes of 
S n+l and the torus, ignoring all the higher harmonics. There will be two contributions to 
the potential of these scalars: one due to the curvature of the S n+1 that we compactify on, 
another due to the kinetic term of the 4-form (3-form) potential in IIB (eleven-dimensional) 
supergravity. Both of these terms will turn out to be exponentials of a linear combination 
of the two scalars, hence the notion of universality classes introduced in section |3] will be 
relevant in this context. 



4.1 Kaluza-Klein reduction 

We want to perform the dimensional reduction in a way that would give an Einstein-Hilbert 
action for gravity, without a scalar-dependent coefficient in front of the three-dimensional 
Ricci scalar. For this purpose, let's start with a general metric 

ds 2 D = ds 2 r + e 2K ds 2 p + e 2L ds 2 q (32) 

where D = p + q + r, ds 2 is a metric on the Lorentzian r- manifold X r , ds 2 and ds 2 are metrics 
on compact Riemannian manifolds Y p and Z q , respectively, and K and L are functions that 
depend only on the coordinates of X r . For example, for the case of the D3-brane, p = 2, 
Y p = T 2 (two worldvolume directions are compactified) , q = 5, Z q = S 5 , r = 3. Defining 

ds 2 = e 2N ds 2 r , (33) 

with 

N = _pK + qL 
r — 2 

one gets 

y/=gER(df D ) = V^FrVFpV^l^l) + e 2 ^R(ds 2 p ) + e 2 ^R{ds 2 q )- 

1 n (35) 
-(p(r + p- 2)d fl Kd v Kg?" + q(r + q - 2)d^Ld u Lg^ + 2pqd fl Kd v Lg^ 



where R(ds 2 ) denotes the Ricci scalar for the metric ds 2 . The first term gives the Einstein- 
Hilbert lagrangian for ds 2 (without a scalar-dependent factor — as we wanted), and the rest 
are the kinetic and potential terms for the scalars K and L. In our applications, R(ds 2 ) and 
R(ds 2 ) will be the Ricci scalars of T p and S q , respectively, i.e. and q(q — 1). 

Using these results, we get the dimensionally reduced form of the D3-brane metric and 
the scalars as 

ds 2 = -r w fhdt 2 + r w ^dr 2 + r 10 fdx 2 
K (r) = ~]ogf (36) 
L {r) = logr + jlog/ 
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where / and h are given by (J2J) with n = 4. The three-dimensional action is given by 



J d 3 x^ (r - \d^d v ^G l3 g» v - (37) 



Ga = ( 11 Z 




S 

where 

v 20 60 

V{$) = -20e~ 4K ~ 12L + 8i% 3 e- 4 *- 20i , (39) 

and Rd3 is the length scale defined in (J7J), with n — 4. The first term in V comes from 
(135)1 . and the second one comes from the dimensional reduction of the _F| term in the IIB 
SUGRA action. We have omitted overall factors from (|37|). and we have introduced notation 
that somewhat obscures dimensional checks: for instance, e L has the dimensions of length. 
In practice, we will simply choose an arbitrary numerical value for Rd 3 for the numerics in 
section l4~2l Another way to think of this is that we choose an arbitrary value for the Planck 
length, which is permissible because our analysis is entirely classical. 

Since the second term in the potential is due to the charge of the D3-brane, we expect 
this term to be of little significance in the large non-extremality limit. By performing a 
linear transformation of the scalars K and L to canonically normalized ones (i.e. those with 
Gij = 5y in equation (|15|)) one of which is aligned with the exponent in the first term of 
the potential, one finds that this term corresponds to a Chamblin-Reall coefficient 7 = 
That this value is in the unstable range (|27|) is consistent with expectations, since the brane 
is expected to be unstable when the effects of the charge are negligible. 

Another consistency check can be done by calculating the relevant exponent a in the 
Chamblin-Reall equation of state by using (|21)l . The equation of state of an uncharged 
brane can be read off from (jHJ) and (jOJ to be S oc VT~( n+l \ The exponent obtained from 
()21j) by setting 7 = y| agrees precisely with this: one gets S oc VT~ 5 in both cases. 

When the brane is near-extremal, R 4 = r 4 sinh 2 a ^> r 4 , so for radii close to the horizon, 
R^> r. Using this in (|3fi|) . we see that Lq is approximately constant near the horizon, i.e. for 
a near-extremal D3-brane the "active" scalar close to the horizon is K. Thus, in order to read 
off the relevant Chamblin-Reall exponent 7 in the near-extremal regime, one sets L = const 
in the action, and rescales K to make it canonically normalized: K = \f\2K. Then, both 

terms in the potential become a multiple of e V3 which shows that 7 = -^=. This value of 7 is 
in the thermodynamically stable range, as expected for a near extremal D3-brane. Equation 
(12 lj) gives the equation of state of the near-extremal D3-brane to be S oc VT 3 , which is the 
equation of state of a 3+1-dimensional CFT, in line with the AdS/CFT correspondence. 

The expected transition from the stability for near-extremal branes to instability at large 
non-extremalities is also confirmed by the numerical analysis of dynamical instabilities, which 
we present in the next section. 



4 The sign of 7 is arbitrary since the field redefinition <f> — > —0 is equivalent to switching the sign of 7 in 
JUJl, HU and 123). 
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The analysis for black M2 and M5-branes is similar. M2-brane versions of (|36|). (|38|1. and 
(JUJ) are 

dsjj = _ r i4 /Mt 2 + r u ^dr 2 + r u fdx 2 

K (r) = -glog/ 
L (r) = logr + ^log/ 




= -42e- 2 ^- 16L + 18Rf I2 e- 2K ~ 28L 



and M5-brane versions are 



d s 2 = ~r 8 fhdt 2 + r 8 ^dr 2 + r 8 fdx 2 



h 

Koir) = --log/ 
o 

Loir) = logr + -log/ (41) 



3 

40 32 
32 40 



V($) = -l2e~ 8K ~ WL + 9 -RW~ 8K - 1GL , 

where h and / are given by (J2J) and i?M2 and i?M5 are given by J7J) with n = 6 for the 
M2-brane and n = 3 for the M2-brane. Once again for large non-extremality the potential 
terms due to the charges can be neglected, and by using canonically normalized scalars one 
can calculate the 7's that are relevant for the terms due to the curvatures of S n+1 . One gets 
7 = -7= for the M2-brane and 7 = ^7= for the M5-brane. As in the D3-brane case, these are 
consistent with the equations of state of uncharged branes: S oc VT~ 7 for the M2-brane and 
S oc VT~ A for the M5-brane. In the near-extremal limits, the "active" scalar is K for the M2 
and M5-branes as well, and one gets the corresponding Chamblin-Reall coefficients as 7 = | 
for the M2-brane and for the M5-brane. The equations of state one gets from these are 
consistent with expectations from AdS/CFT: S oc VT 2 for the M2-brane and S oc VT 5 for 
the M5-brane. 



4.2 Numerical stability analysis for non-dilatonic branes 

We will describe the stability analysis for the D3-brane in some detail, and just show results 
for the M2-brane and the M5-brane. 
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In order to investigate the dynamical stability of the three-dimensional background (J36j) . 
obtained by Kaluza-Klein reduction from the D3-brane, we introduce the perturbations 



ds 2 = _ e 2A(r, X ) (rlOf^ + ^(r,*) j ^0^2 + ^ f ^2 



h 



A(r,x) = A cos kxA\{r) 
B(r,x) = A cos kxB\ (r) 
K(r,x) = K (r) + A cos kxK^r) 
L(r, x) = Lq{t) + A cos kxL\{r) , 



(42) 



and work to linear order in A. K and L are the scalars introduced in (|37|). 

There are four nontrivial Einstein equations (those involving R\, R r r , R x x and R r x ), and two 
scalar equations of motion. Using the R r x equation, one can solve for B\ = B\(Ai, A[, K\,L\). 
The R\ equation and both of the scalar equations of motion involve Bi, but not its derivatives. 
Plugging in B\ as found from the R r x equation into the R\ equation and the scalar equations 
of motion, we end up with three equations involving A\, K\ and L\ (it is possible to do a 
consistency check by arriving at the same equations using the R x and R r r equations). Defining 
R 4 = rQsinh 2 o;, the equations read 



At(r) 



=r \16 r 2 r 4 (5 r 8 + 10 r 4 R A + 3 R 8 - 2 R A r 4 ) 

+ k 2 (r 4 + R A ) 2 (5r 8 - i? 4 r 4 + 3r 4 (r a - r 4 )) \ 
-32r 3 r 4 (5r 8 + 10 r A R A + 3 R 8 - 2 R A r A ) K x (r) 

- 160 r 3 (2 r 4 R A r A + R 8 r A + r 8 (2 R A + 3 r 4 ) ) L x (r) (43) 

- 2 ^5 r 16 + R 8 r 8 - 5 r 12 (i? 4 + r 4 ) + r 8 (20 i? 4 r 4 + 6 r 8 ) 

+ r 4 (S^^-Si^ro 8 ) J A\{r) 

- r (r 4 + # 4 ) (r 4 - r 4 ) (5 r 8 - i? 4 r 4 + 3 r 4 (R A - r 4 )) A*(r) 
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=8 r 3 r 4 (lO (r 4 + R 4 ) 2 - 4 i? 4 (i? 4 + r 4 ) ) Ax (r) 
+ ^-32 r 3 r 4 (5 r 8 + 10 r 4 i? 4 + 3 # 8 - 2 i? 4 r 4 ) 

+ 3A; 2 r(r 4 + J R 4 ) 2 (5 r 8 - R 4 r 4 + 3 r 4 (i? 4 - r 4 ) ) W (r) 
+ 5 (V r (r 4 + i? 4 ) 2 (5 r 8 - i? 4 r 4 + 3 r 4 (# 4 - r 4 )) 

-32r 3 (2r 4 R 4 r 4 + R 8 r 4 + r 8 (2 R 4 + 3 r 4 )) ^ L x (r) 

+ (r 4 -r 4 ) (15 r 12 + 3 R 8 r 4 + 5 r 8 (lO i? 4 + 3 r 4 ) + r 4 (l5 R 8 - 2 i? 4 r 4 )) A'^r) 
-3 (r 4 + i? 4 ) (5r 4 -r 4 ) (5 r 8 - R 4 r 4 + 3 r 4 (# 4 -r 4 )) tf((r) 
-5 (r 4 + i? 4 ) (5r 4 -r 4 ) (5 r 8 - i? 4 r 4 + 3 r 4 (# 4 -r 4 )) Li(r) 

- 3 (r 4 + R 4 ) (r 5 - rr 4 ) (5r 8 - R 4 r 4 + 3r 4 (i? 4 - r 4 )) A^(r) 

- 5 (r 4 + i? 4 ) (r 5 - rr 4 ) (5r 8 - R 4 r 4 + 3r 4 (i? 4 - r 4 )) L'/(r) 

(44) 

=8 r 3 r 4 ^6 (r 4 + R 4 ) 2 - 4 i? 4 (i? 4 + r 4 ) ) A a (r) 
+ ^ - 32 r 3 r 4 (3 r 8 + 6 r 4 i? 4 + i? 8 - 2 i? 4 r 4 ) 

+ fc 2 r(r 4 + i? 4 ) 2 (5r 8 -i? 4 r 4 + 3r 4 (r 4 - r 4 )) ^ #i(r) 
+ f3fc 2 r(r 4 + iZ*) 2 (5 ^-^^ + 3^ (# 4 -r 4 )) 

- 32 r 3 (3 i? 8 r 4 + r 8 (6 R 4 + 9 r 4 ) + r 4 (-4 R 8 + 2 R 4 r 4 ) ) L x (r) (45) 

+ (r 4 - r 4 ) (9 r 12 + _R 8 r 4 + r 8 (34 A 4 + 9 r 4 ) + 5 r 4 (i? 8 - 2 _R 4 r 4 ) ) A[ (r) 

- (r 4 + R 4 ) (5 r 4 - r 4 ) (5 r 8 - i? 4 r 4 + 3 r 4 (i? 4 - r 4 ) ) AT( (r ) 

- 3 (r 4 + R 4 ) (5 r 4 - r 4 ) (5 r 8 - i? 4 r 4 + 3 r 4 (# 4 - r 4 ) ) Li (r) 

- (r 4 + i? 4 ) (r 5 - rr 4 ) (5r 8 - i? 4 r 4 + 3r 4 (r 4 - r 4 )) iTf(r) 

- 3 (r 4 + R 4 ) (r 5 - r r 4 ) (5 r 8 - i? 4 r 4 + 3 r 4 (R 4 - r 4 )) L'[{r) . 

Solving these linear ODE's numerically, one extracts using a shooting algorithm the 
value of k where a static perturbation exists, for any specified non-extremal mass. The 
results of this numerical investigation are summarized in figure 121 For large non-extremality, 
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Figure 2: Instabilities of black D3, M2, and M5-branes. The dots represent threshold insta- 
bilities of the black branes we found. The vertical lines indicate the thermodynamic stability 
limits, which were found by JEJ, with R = 1 for the D3-brane, R = (1/3) 1 / 6 for the M2- 
brane, and R = (2/3) 1 ^ 3 for the M5-brane. The curves indicate wavenumbers of threshold 
instabilities of uncharged black branes, with tq being the Schwarzchild radius. 
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the threshold wavenumbers fit nicely onto the curve that gives those of an uncharged 3- 
brane. (This curve was obtained by reading off the wavenumber of the threshold instability 
of an uncharged 3-brane from the plots in |2], and doing an appropriate scaling to get 
the values for an arbitrary ro.) As one approaches the thermodynamic stability limit, the 
wavenumbers approach to zero. In fact, by zooming into this region, we were able to show 
that the wavenumbers go to zero approximately with a power-law behavior, with the critical 
exponent being close to 1/2: that is, k ~ V r o — r o where Tq is the critical value where 
the instability disappears. Because M is a non-singular function of ro, one could also write 
k ~ y/M — M* . A naive way of understanding this is to think of some effective theory 
describing the perturbations in terms of a bosonic field with m 2 < 0. Then the criterion for 
the perturbations to be static is uj 2 = k 2 + m 2 = 0, where k is the wave-number of a static 
perturbation. The resulting equation, k = \f—m 2 , predicts the observed scaling, provided 
m 2 ~ Tq — r . Such analytic behavior is a standard assumption: our argument basically 
amounts to a naive application of Landau theory. 5 

The numerical stability analyses for the M2-brane and M5-brane are similar to the one for 
the D3-brane. The results are summarized in figure 121 Once again, near the thermodynamic 
stability limit, the wavenumbers go approximately as y/r — rj. 

5 Conclusions 

The stability of charged p-brane solutions was first attacked in [2] . Although the methods 
presented there are in principle practicable for any solution, they require considerable com- 
putational fortitude to apply. We have argued that a simple notion of universality classes 
is fairly broadly applicable to the question of brane stability: one commonly finds a ther- 
modynamic relation S oc VT a for charged p-brane solutions in some limiting regime of 
parameters, and when one does, the stability properties depend only on a — stability per- 
taining when a > 0. Given our treatment of the problem via a Kaluza-Klein reduction / 
truncation to a 2+1- dimensional action involving scalars and gravity, the dispersion relation 
of unstable modes would also appear to be the same for any solution in the universality class 
labeled by a given a. 

We have further shown, in section 01 that when one approaches a boundary of stability, 
the critical wavelength separating stable from unstable modes diverges, and that it does so 
for the cases studied with a critical exponent that could be guessed on the grounds of a 
naive effective field theory argument. It would be interesting to see if other exponents arise 
from non-generic situations, for instance a case where the temperature T depends on mass 
as T max — T ~ \M — M*\P for some [3^2. It is possible that spinning brane solutions 
provide a venue for more intricate thermodynamics [2U I2H] , but one would encounter there 

5 It is actually a little too naive, because it predicts an incorrect dispersion relation. In EEL ^ was found 
that uj — > as k — > 0, whereas using ui 2 = k 2 + m 2 would suggest finite u> at k = 0. Possibly an improved 
understanding could be based on hydrodynamic considerations. 
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the complication of chemical potentials for the various angular momenta. 

For the p-branes of type II string theory and M-theory, away from extremality but without 
angular momenta, thermodynamic stability pertains up to an upper mass limit given in 
equation (jllj) . For the D3-brane, the M2-brane, and the M5-brane, our numerical analysis 
supports the claim [21111 E] that dynamical and thermodynamic stability coincide. Extending 
the analysis to other p-branes should not be too difficult. 

In conclusion, it seems that the notion of universality classes that we introduced in sec- 
tional together with the behavior at a boundary of stability explored in section 01 represent 
a fairly comprehensive description of the qualitative features of the Gregory-Laflamme insta- 
bility in linearized perturbation theory. Of course, it would be desirable to go beyond clas- 
sical perturbation theory in understanding universality classes of behaviors for non-uniform 
horizons. Some of the methods of this paper may prove useful in that broader context as 
well. 
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